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1. INTRODUCTION 

Type | diabetes is an inveterate disease in which the pancreas does not produce sufficient insulin to 
control glucose levels in the blood. In the shortage of insulin, the admittance of glucose into skeletal, cardiac, 
smooth muscle and other tissues is reduced. When insulin is missing for a longer period of time, the muscle 
and tissue cells will start using fat as energy source, in place of glucose from the blood stream. 
The disturbance of meal eating presents great challenges to automatic blood glucose control. There are many 
factors must be measured to regulate blood glucose such as. Fats and proteins cause delays in absorption of 
glucose from carbohydrates eaten at the same time, physical exercise affects the blood glucose regulation by 
reducing the requirement for insulin, and furthermore, stress and physical illness are involved in the blood 
glucose regulation process. Thus, designing accurate predictive models are hard to achieve. In 2011, Dimitri 
Boiroux et al, applied a robust feedforward-feedback control strategy to people with type | diabetes. The feed 
forward controller consists of a bolus calculator which compensates the disturbance coming from meals. 
The feedback controller is based on a linearized description of the model describing the patient [1]. 

In 2008, Gianni Marchetti et al concerned with the improvement of new feed forward-feedback 
control techniques for actual glucose control and type Idiabetes. Enhanced post-meal responses can be 
achieved by a pre-prandial snack or bolus, or by reducing the glucose set-point prior to the meal [2]. In 2017, 
Mojgan Esna-Ashari et al proposed a technique for blood glucose level regulation in type I diabetes. 
The control approach is based on non-linear model predictive control. The purpose of the proposed controller 
optimized with genetic algorithms is to measure blood glucose level every time and estimate it for the next 
time interval [3]. In 2017, Lukas Ortmann et al, introduced a Model Predictive Controller that takes 
the periodic Insulin sensitivity into account, so as to improve blood glucose control. The upcoming effect of 
the Insulin sensitivity is estimated by a machine learning technique, namely, a customized Gaussian Process 
(GP), based on historical training data [4]. In 2017, Chiara Toffanin, proposed an adaptive Model Predictive 
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Control approach to account for it and test it in silico. Run-to-Run method adapts the subcutaneous basal 
insulin delivery throughout the night and the carbohydrate-to-insulin ratio throughout the day, according to 
some performance indices measured from subcutaneous continuous glucose sensor data [5]. In 2017, Berno J. 
E. Misgeld et al presented a novel method towards real-time estimation of insulin sensitivity in type-1 
diabetes. The predictive of insulin sensitivity is investigated while a carbohydrate uptake disturbance takes 
place simultaneously [6]. In construct, in this paper, the Robust Model Predictive Controller (RMPC) design 
approach will be employed for the regulation of constrained systems. RMPC algorithms are able to 
coordinate hard constraints on state outputs and inputs. Stability and robustness are considered properties of 
RMPC controllers, as they are essential to hardware implementation. A whole theory on robust stability of 
RMPC has been applied for linear and nonlinear systems. 


2. PROBLEM DESCRIPTION 

The controllers are built in both traditional and RMPC. However, it holds only for linear scheme of 
the system. While dealing with the non-linear system operations, the system should be linearised. 
This linearised form of the non-linear system could be properly control using traditional controllers such as 
(i.e. Model Predictive Controller). In an Artificial Pancreas regulation, continuous subcutaneous insulin 
infusion common with subcutaneous glucose sensing has been shown to be effective for feed forward 
channel “out-of meal” periods, preservation blood glucose in the regular level in the post absorptive case. 
However, the feedback channel has not been so successful in treating insulin desires at meals state [7, 8]. 
In spite of developments in sensor and pump equipment, the Artificial Pancreas must coordinate with 
the delays and inaccuracies in glucose sensing. This is very hard when a system disturbance, such as, a meal, 
happens and triggers a fast glucose increase that is considerably faster than the time required for insulin 
absorption. The stability problem of constraints plants in closed-loop with MPC controllers is shown in 
Figure 1. 
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Figure 1. Artificial pancreas techniques 





3. METHODOLOGY 

Model predictive control is a method of control in which the present control signal is obtained by 
solving, at each sampling time, a finite constraints open-loop optimal control problem, using the recent state 
of the system as the initial state, the optimization produces an optimal control signal sequence and the first 
control signal in this sequence is applied to the system. MPC is able to cope with hard constraints on control 
signals and states as shown in Figure 2. Therefore, been applied in medical application where approval of 
constraints is mostly very significant since efficiency requests operating points on or close to the limit of 
the set of acceptable states and control signals. 


Onniirhance 
Actual Bödy 


Blood glucgse 


f level 
Estimator z Constraints = 


system 





mor 
Bounda 


Figure 2. MPC for constraints system 
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4. MODELING AND SIMULATION 

In this paper, the MPC design will be employed for the regulation of constrained systems. A robust 
stability of MPC has been avoid violations of input and output constraints drive certain output variables to 
their optimal set points, while keeping other outputs within specified value. Avoid redundant movement of 
the input variables. Estimated values of output variables are predicted using a linear quadratic regulator of 
the plant and current measurements. The control signal calculations are according to both upcoming 
predictions and recent measurements. The control signalsu(k), at the kt” sampling time are calculated so that 
they minimize a performance index J. the constraints of optimal control system as (1): 


ee <u(k) < u(k) max (1) 
Y(k)min < y(k) = V(kK) max 


The Insulin kinetic equations are (2): 


i =—[-I(t) + kS) 
ž = [-X®) +10) (2) 
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Where, I- insulin plasma, X- insulin interstitial, S- insulin sensitivity and U- control signal 
While, the glucose kinetic equations are (3) [9]: 
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where, G- glucose plasma, and Y- subcutaneous compartment. 
According to medical experience Ty, ki, Tm, Ti, Ty, and Tg can be considered constant. According to 


insulin and glucose kinetic equation the state space representation for both insulin and glucose equations are 
(4-5): 
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According to medical analysis the hepatic balance equations was designed to reduce the effective of 
meal (disturbance) (6) 


e, if 0<e,<0.88 


0.88 if e, = 0.88 
da = | 
0 otherwise 


qg = 0.06 G — 0.25 


af=fqg if qg 20.25 
0 otherwise 
dc = 0.25E, + qg (6) 
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en, = 0.87 + tanh0.0045(G — 175) 
Cr = 0.11e,, 


From the constraints system output, the assumed input and output are given as known parameter to 
the constraints system identification toolbox. It provides the best appropriate transfer function by using 
the adaptive tuning method to get the values of poles and zeroes and linearized the constraints system which 
is shown in Figure 3. The achieved system model is used to linearise the constraints system and the model 
predictive is designed to simulate the system as presented in Figure 4. 
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Figure 3. Measured and simulated system output 
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Figure 4. Input state (insulin) to output state (blood glucose) for plant 


After finding the traditional (PID) controller, the plant is obviously checked and estimates 
the performance. However, the settling time taken is infinite and the process is not stable for the estimated 
output. While, when the RMPC is connected vice PID controller. The RMPC that have been designed 
the controlled variables and model output based on MPC controller toolbox. The inherent designed RMPC 
controller had model which is shown in Figure 5. RMPC can increase system speed and estimation accuracy. 
This technique gives a fairly proper estimate of the speed and is robust to any parameter variation. 
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Figure 5. RMPC output for insulin and blood glucose states 
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5. CONCLUSIONS 

In this paper, we described RMPC based on a linearized form of constraints system. The state is 
increased with an inserted meal (disturbance) to ensure a state control of the blood glucose. We use MPC 
techniques to calculate insulin quantity in the cases where the insulin sensitivity reductions and in the case 
where the any meal (disturbance) is taken. The simulation results presented in this paper show that 
the designed RMPC is more effective, it has faster settling times. The uncertain parameter in input and its 
effect on the different factors of the real-time system are also measured. The designed RMPC not only 
considered of the unexpected disturbances in system, however also recover the certain parameters to 
the desired values, hence presenting the robustness in performance. 
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